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Abstract When a gravitating object moves across a 
given mass distribution, it creates an overdense wake 
behind it. Here, we performed an analytical study of 
the structure of the flow far from object when the flow is 
isentropic and the object moves subsonically within it. 
We show that the dynamical friction force is the main 
drag force on the object and by using a perturbation 
theory, we obtain the density, velocity and pressure of 
the perturbed flow far from the mass. We derive the ex- 
pression of the dynamical friction force in an isentropic 
flow and show its dependence on the Mach number of 
the flow and on the adiabatic index. We find that the 
dynamical friction force becomes lower as the adiabatic 
index increases. We show analytically that the wakes 
are less dense in our isentropic case in comparison to 
the isothermal ones. 

Keywords hydrodynamics - ISM: general - galaxies: 
kinematics and dynamics - stars: kinematics 



1 Introduction 

The motion of a gravitating mass accross a given dis- 
tribution of matter creates an overdense wake behind 
it. In turn, this overdense region interacts gravita- 
tionally with the moving object causing it to lose mo- 
mentum and kinetic energy. The ambient medium 
gains the momentum and the kinetic energy lost by 
the moving mass. The drag force exerted on the 
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moving obj ect is also known as the dynamical fric- 
tion force. IChandrasekharl (J19431 ) studied the effects 



of the dynamical friction force in the collisionless sys- 
tems analytically and his results has been applied to 
many astrophysical systems. The effects of the dy- 
namical friction force have been investigated in the 
context of the evolution of massive you ng stellar clus 



ters near to the Ga lactic center (e.g.. iGerhardl 12001 



Kim fe Morris 20031). the dynamical evolution of s tellar 



clusters (e.g.. IPortegies Zwart fe McMillan! 120021) . and 
the migration of planets in a protopl anetary disc made 
of collisionless planetesima ls (e.g., iDel Popolol 12003 : 
Muto. Takeuchi fe Idall201ll ) . N-body simulations have 



also shown that the dynamical friction force based on 
the Chandrasekhar's formula c an be extend ed to the 
many-body interacting systems ( Inouell2011 ). 

Dynamical friction effects are also present when 
a gravitating mass moves accross a collisional sys- 
tem, i.e. a gaseous medium. The work done by 
the dynamical friction for ce may act a s a he ating 
mechanism of the gas. lEl-Zant et al.l (|2004T) 



Faltenbacher et al 



and 



(|2005f ) have shown that the in- 
tergalactic gas in the central regions of the galaxy 
clusters can be heated by the motions of the galax- 
ies. Other authors have investigated the role of the 
gas drag force in the context of pla net migration in 



Ida fe Lin 1996 



gaseous protoplanetary di scs (e.g , 
Haghighipour fe Boss! 12003 ). the dynamics of the pro 
tostars and the sta rs in the stellar clusters prior to the 
gas removal (e.g.. lEscala. Larson. Coppi fe Mardones 



gas removal (e.g. , Hiscaia. Par son, cop 
20031 ISaivadpour. Deiss fe Kegell Il997l ) 



and the mo- 
tion of bla ck holes on circular orbits in a gas rich galac- 
tic nuclei JKim fe Kimll2009l: lKimll2010[ ). The way gas 
drag affects the collisional behaviour of planetesimals in 
a gaseous di sc has also been investigated by several au- 
thors (e.g., Kobavashi. Tanaka. Krivov fe Inabal bOlOt 



Ormel fc Klahrll2010D . 



From a more fundamental point of view, a num- 
ber of authors studied the nature of the drag force 
when the perturber moves supersonically in the gaseous 
medium ( Dokuchaevj 1964t iRuderman fc Spiegell Il971[ 
Rephaeli fc Salpeterl 1 1980). In the s teady-state sub- 



sonic case, 



Rephaeli fc Salpeterl (119801) co ncluded that 



the drag force is absent. lOstrikerl ( 1999 ) studied the 
time-dependent linear theory of the dynamical friction 
force and found the value of dynamical friction is not 
zero in contra ry to the previous studies. Simu l ations 
performed by ISanchez-Salcedo fc Brandenburg! (119991) 



confirmed the analytical results of lOstrikerl (|1999l ). al- 
though other simulations find a different result (e.g., 
Ruffert 1996). 

More recently, iLee fc Stahlen |201ll ) (hereafter LS) 
considered the role of the dynamical friction in the sub- 
sonic case. They argued that because time does not 
appear explicitly in Ostriker's solution, a steady-state 
analysis can be applied in the subsonic case. In fact, the 
dynamical friction force should reach a steady state if 
the object travels for a sufficiently long time through a 
very extended medium. In all simulations of this prob- 
lem, a steady-state value seems to be approached. It is 
true that, in linear perturbation theory, there is no dy- 
nami cal friction force in the subsonic regime (jOstriker 
1999). There is also no mass accretion onto the ob- 
ject in linear theory, as was shown in LS. Of course, 
neither the force nor the accretion rate is truly zero. 
Both arc finite if one extends the perturbation analysis 
to second order. Also, LS assumed the object's radius 
is much smaller than the accretion radius as it is truly 
valid for real astronomical objects and under this as- 
sumption the relative density enhancement in the wake 
is not small and so they could extend their analysis to 
the nonlinear case. 

Ostrikcr (1999) attempted to circumvent this zero- 
drag problem by assuming that the gravitational force 
from the moving body suddenly switches on at a cer- 
tain time. There is no physical basis for this assump- 
tion. Moreover, Ostriker calculated the force by in- 
tegrating over linear density perturbations in the gas. 
But, under her assumption that the physical size of the 
object is much smaller than the accretion radius, the 
density perturbations are nonlinear close to the object. 
LS approach uses perturbation theory where it is indeed 
valid, far from the object. We also accounted for that 
portion of the force arising from direct accretion of mo- 
mentum, a contribution neglected by Ostriker. On the 
other hand, not all the numerical simulations have con- 
firmed Ostriker's result, in particular those which as- 
sumed a small accretor obtained a very different force 
(e.g., Ruffert 1996). But LS results showed that the 
friction force is given by MV (where M is the mass 



accretion rate onto the object/perturber and V is its 
velocity) as long as the perturber is moving subsoni- 
cally. 

In the present work, we extend the approach of LS to 
the isentropic case. We derive the steady-state analyt- 
ical solutions for the dynamical friction force when the 
ambient medium is isentropic. Our results show that 
the dynamical friction force remains equal to MV sim- 
ilar to the isothermal case but in our case M is the mass 
accretion rate when the gas is isentropic. We focus on 
the region far from the object and we also assume that 
the object's radius is much smaller than the accretion 
radius. The basic equations are derived in § 2. The 
perturbed equations and the boundary conditions are 
introduced in § 3 as well as the mass accretion rate and 
its relation to the stream function are obtained in this 
section. In § 4, wc solve the dynamical equations for 
the first-order perturbed equations. In § 5, we solve 
the second-order equations and the dynamical friction 
force is calculated in the isentropic regime. Our results 
arc summarized in § 6. 



2 basic equations 

The basic equations of our problem in the steady state 
can be written as 



V • (pU) = 0, 



p(U-V)U = -VP-p 



(U.V)£?-4(U.V> 

p 



GM 



o, 



(i) 



(2) 



(3) 



where p, U, P, E are the density, the velocity, the 
pressure, and energy, respectively. The energy E is 
defined as E = 1/(7— l)(P/p), where 7 is the adiabatic 
index. Here, R is the radial distance in the spherical 
coordinates (R, 8, (p) whose origin is anchored on the 
gravitating mass M . We also assume that the gas is 
isentropic with 7 > 1. Wc assume that very far from the 
mass which travels in a straight line in the medium with 
velocity V, that density is uniform, i.e. lrnifl-^xj p = po- 
It is more convenient to introduce a set of dimcnsionless 
variables which are defined as r = R/R s , u = U/c s , 
g = p/po and p = P/Pq, where c s is the sound speed 
and R s is the sonic radius which is defined as R s = 
GM/c^. Equation Q is the energy equation and we 
also introduce the dimensionless variable e = E / Eq = 
p/g. Using these new variables, equations (H}-([3|) can 
be rewritten as 



V • (gu) = 0, 



(4) 
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g(u.V)u = -Vp-4, (5) 



(u.V)^ l -^(u.'V)g = 0. (6) 

It is clear that, in equations (HJ)-©, the radial deriva- 
tives are calculated with respect to the non-dimensional 
radial distance r. We also define j3 = V/c s , where V 
is the velocity of the flow at the large distance from 
the object. Thus, properties of the flow at the dis- 
tances far from the gravitating mass M are u = /3e z , 
g = 1, p = 1, when r>l. Note that in our calcu- 
lations it is assumed that the system is axisymmetric 
(d/ 'dip = 0). It is advantageous to introduce the non- 
dimensional stream function -0(r, 9) which is defined as 
being 



gu = V x ( 



* 



rsin9 



(7) 



then the non-dimensional continuity equation (j4]) is au- 
tomatically satisfied. Thus, the components of the ve- 
locity can be written in terms of the stream function as 



1 



dip 



gr 2 sin 6 89 



and 



u e = 



-1 dip 

grsiniO dr 



(8) 



(9) 



Having the above equations for the velocity in terms 
of stream function ip, we are left only with the momen- 
tum equation ([5]) and the energy equation ([5]). Thus, 
we have 



du T 



dr 



ug du 
r ~d~6 



dp 
dr 



Q_ 

,•2' 



(10) 



g u T 



dug 

Or 



Ug dug U x Ug 

r d9 r 



Mr 7^ + - -^ I «r^ ^ 



dr 



ug de 
Vdd 



1 dp 
rd9' 

c)g 



dr 



ug dg 
r d9 



(12) 



The above equations constitute our main dimensionlcss 
equations. In the next section, we perturb them using 
perturbation expansions. 



3 perturbation expansions 

The stream function and the density are expanded as 
in LS, i.e. 



4(r,0) - / 2 (% 2 +/i(9)r+/„(«)+/-i(«)r- 



(13) 



g{r, 0) = l + g-x{e)r- x + . 9 - 2 (#K 2 + ff _ 3 (e)r~ 3 + ■ 



(14) 



and in our isentropic case, we need to expand the pres- 
sure as well 



p{r, 6) = l + h-^ey- 1 + h^ 2 (8)r- 2 + fc_ 8 (0)r 



-3 



(15) 



Considering the fact that velocity is constant far 
from the mass, fziQ) can be obtained as /2(d) = 
(l/2)/3sin 2 6. Since we set tp(r,Tr) = 0, then fi(ir) = 0, 
for i = 1, 0, — 1, —2, ect. The behaviour of the velocity 
gives us more boundary conditions as /|(7r) = //(0) = 
fori= 1,0, -1,-2, ect., and /j(0) = for i = 1,-1,-2, 
ect. (for further details see LS). Substituting the above 
series expansions into the main equations, we obtain a 
set of ordinary differential equations. 

We can now derive a relation between the accretion 
rate and the coefficient of the density expansion. From 
the continuity equation (fTJ), the mass accretion rate is 
written as 



M= -2tt / pu R R 2 sin8d6. 
'0 



(16) 



Substituting the expansion for the stream function 
into equation ([5]) for the radial velocity u r , the non- 
dimensional form of M in units of 2irpoc s R 2 becomes 
(see LS for details) 



(11) M = / o (0). 



(17) 



From equation ([TT]) , we can see that the mass accretion 
rate depends just on the coefficient /o at 9 = which is 
a term that is independent of the radius r in expansion 

disD. 



4 First-order equations 

We can now substitute the expansions (flU|). (fLf]) . and 
(fl~5|) into the equations (TTU)) and <jTTJ) and the energy 
equation (fj"2"j) . Then, we match the coefficients of each 
power of r. We find that the highest power of r is r _1 , 



but all the corresponding coefficients on both sides of 
the equations vanish identically. So, we consider the 
next power of r (i.e., r~ 2 ) and the first-order equations 
arc obtained. From the radial component of the mo- 
mentum equation, we obtain 

-Pf"-/3fi+P 2 sin6cos6g'_i+P 2 cos 2 6g-i-h-i+l = Oj 

(18) 

and the 9— component of the momentum equation gives 
P 2 sin 2 dg'_! + P 2 sin 9 cos dg- X -h'_ 1 = 0, (19) 

and finally, the energy equation (fT2j) yields 



7 sin 



Considering equation (|18p and the above equations 
for g_i and h—i, a differential equation for f\ is ob- 
tained 



2 °"'_ 1 - sin 2 9h'_ 1 +sm9cos9(-fg-i - h-i) = 0. 



(20) 



We solve the above differential equations analytically 
for /i, g_i and h—i. To the first-order, the solutions 
determine properties of the flow around the moving 
object. We can integrate equation (f2"0")) simply by re- 
arranging its terms as 

sin6>cos6>(7<7_i-/i_i) + sin 2 6»— (jg_i~h-i) = 0. (21) 

d9 

Introducing W — 7.9-1 — h—i, this equation becomes 
dW/d9 = ~Wcot9 and its solution is W = W /sin9 
where Wo is an arbitrary constant to be determined 
from the boundary conditions. Thus, 



7.9-1 — h—i = Wo/ sin ( 



(22) 



From this equation, one can simply calculate h'_ 1 in 
terms of g'_ 1 which enables us to re- write equation ([20)) 
as follows 

{P 2 sin 2 9~ 1 )g'_ 1 +f3 2 sin 9 cos 6g-i- W cos 9/ sin 2 9 = 0.| 

(23) 

This equation is integrable and gives us 
W C 



,„,, 1 lC ( 7 -/? 2 ) 
h + n = ~b 



P /3 (7 - /3 2 si 



2 oi-n^ f)\o 



(26) 



The general solution of the above equation is written 

as 



/i = \--n\Jl- P 2 sin 2 9 + Acos9 + B sin 9. (27) 
P P 

The boundary conditions help us to determine the con- 
stants A, B and C. Since we have fi(w) = and 
/i(0)=0, then 



1 - Ap = CV7, 
and 

1 + AI3 = CV7- 



(28) 



(29) 



Therefore, the constants become C = 1/^/7 and A = 
0. Similarly, the boundary conditions f[(ir) = and 
/i(0) = give B = 0. So, the coefficient fi{9) can be 
written as 



and also we have 

1/V7 



g~m = 



h-i{0) 



V7 - P 2 sin 2 9 



V7 



(30) 



(31) 



(32) 



9-i(d) 



7sin6» y J 1 _ ^2 sm 2 



(24) 



where C is a constant of integration. The above equa- 
tion is valid for the any value of 9 including 9 = tr, but 
equation (|24[) goes to infinity when we set 9 = tt unless 
Wo = 0. Substituting expression of g_\ into equation 
|), we obtain /i_i as 



yj 1 - P 2 sin 2 

Up to now, we have derived the first-order density, pres- 
sure, and stream function perturbations analytically for 
an isentropic flow. For the isothermal case (i.e., 7 = 1), 
our first-order solutions reduce to the first-order solu- 
tions of LS. 



5 second-order equations 

Now, we equate the coefficients of the next power of r in 
the components of the momentum equation and the en- 
ergy equation in order to obtain the second-order equa- 
tions. From the radial component of the momentum 
equation, we obtain 



h-x(9) = 



Ci 



V7 - /3 2 sin 2 I 



(25) 



-Pfo - P cot df + P 2 cos 9 sin 6g'_ 2 + 



2P 2 cos 2 



- 2h-2 = Ax + A 2 + A 3 + A 4 , 



(33) 
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g 2 (7i) = -0.88. 

a - ' p = 0.5 ■ 

,1.2 




°- 4 e/, - 6 



Fig. 1 Profiles of fo(0) and <?-2(#) for different values of 
7 with g-a(ir) = -0.88. 



-4 

-6 

-8 

1.2 

0.6 

^0.0 
1 

-0.6 
-1.2 



D) 



-1.8 



1 I ' I 


^-—z^^^^^ 


-Jj2____^^ 


^^s^yr 


" 1.2___^-- " 


s*^ / 


- ]-*~^^^ 


/ 


: /^ 


9>) = 1.12 - 


- 511^-^ 


p = 0.5 - 


- ' 1 ' 1 


i i i i i _ 


^^^Z^ ^ 


: i=}___^^ 


/// 


i - 2 — — — "^ 


i,i, 


r i.4___—-^ 


5/l___^-^ 


i , i 



0.0 0.2 



M 8fe M 



0.8 1.0 



Fig. 3 Profiles of /o(<9) and g-i(&) for different values of 
7 with 17-2(1") = 1.12. 
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Fig. 4 The density contours correspond to 7 = 5/3 (solid) 
Fig. 2 Profiles of / o (0) and fl _ 2 (0) for different values of and 7 = x (dashed) for the cases with ^ = 5 and g _ %{ ^ = 

7 with 17-2(71") = 0.12. i^2 



where all the terms on the righthand side depend on 
the first order variables as 



1 „:„2 n 



sin 3 8 sin 2 8 sin 2 i 



A 2 = Phg-i - ipfg-i cot - /3/i^x cot 8- 



0fi3 , -i+Pfi'g-i, 



A3 = -3.9-1, 



A4 = 2g-ih-i. 



(34) 



(35) 



(36) 



(37) 



After substituting the expressions of g~i, ft— 1, and /1 
into the right-hand side of equation ([42]) and consid- 
ering equation ([43]) . we can integrate equation ([42]) as 
follows 



/3 2 sin 2 0(ft_ 2 - 7.9-2) = i 



/3( 7 - /3 2 sin 2 



+ K, (44) 



where _KT is a constant of integration. This equation is 
valid over the entire range of 8 including 8 = 0. Sub- 
stituting 8 = into equation (|44[) , one can obtain 



Also, there is another equation from the 8— component 
of equation of motion, 



1 7(7-1) 
2 /? 7 

and so 



K = 0, 



K 



~Pfo ~ P sin 6 9-2 + h '-2 - 2 P cos 8 sin ( 

= 6!+6 2 +B 3 , 

where 

„ _ f 2 ^0 A/{ 



S 2 = /3/iff_i cot0 + /3/is_i + 2/3/^1!, 



7-1 
2/3 ' 



£ 3 = -2 5 _ 1 ft'_ 1 . 



(38) 



(39) 



(40) 



(41) 



Now, equation ([4"4"]l becomes 

7-1 1 



ft-2 - 75-2 



2/3 ( 7 -/3 2 sin 2 0) 



(45) 



(46) 



(47) 



Substituting the analytical expressions of /1 and g-i 
into the righthand side of equations ([33]) and ([38]) and 
using equation (|47|). the following coupled differential 
equations for f$ and g_ 2 are obtained 

-Pfo - P cot 6>/^ + /3 2 cos 6» sin 8g'_ 2 + (2/3 2 cos 2 (9 - 27)g_ 2 | 



Finally the last equation is obtained from the energy 
equation, 



1 +±(vl- J_ 

nr. \ 7-1 PFZ 



/3 D V7 V D ^D Ji + VD' 



(48) 



and 



/3 sin 2 0ft'_ 2 - 7/3 sin 2 %1 2 + 2/3 sin cos 0ft_ 2 - 2 7 /3 sin f| _ /3/ , + (/y _ ^ gin 2 ^ _ 2/3 2 cog gin ^ = ^2 gin CQS ^ 
x cos8g-2 = (7 — l)/3sin#cos#<7_ift_i + 7/3 sin 2 ^g^ft-il 



K [3D 2 



2 
D 2 



-/3sin 2 eg^hU - /ift'-i - ft-i/{ + 7.A.9-1 + 7/i.9-i- 

(42) 

The right-hand side of the equations ([33]) , (f38|) and 
(|42p are determined analytically because we have al- 
ready found the coefficients / 1? g_! and ft_i analyti- 
cally. 

One can rewrite the left-hand side of equation ([42]) 
as. 



(3sin 2 d— (ft_ 2 -7.9- 2 ) + 2/3sin6»cos6>(ft_ 2 - jg-2) = 
do 



^d\ lD 7 (^ + /d) 2 



), (49) 



where D = 7 — /3 2 sin 2 0. 

Before solving the second-order equations numeri- 
cally, we obtain a relation between the mass accretion 
rate and the second-order coefficient of the density ex- 
pansion. The left-hand side of the second-order equa- 
tion (1491) is re-written as 



(50) 



-Pfo + (7 - P 2 sin2 %-2 ~ 2 P 2 cos e sin 89-2 



d 



Te (-Pfo + Dg_ 2 ). 



^[/3sm 2 0(ft_ 2 - 75 _ 2 )]. 



(43) 



Having the above equation for the left-hand side of 
equation ([50]) and the fact that D is an even function of 
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9, we can integrate both side of the equation (|49[) from 
9 = -k to which yields, 



(-/3f +Dg_ 2 ) e 



(-13 f + Dg_ 2 ) 



9=0 ' 



(51) 



where the right hand side of the equation vanishes and 
since fo(n) = and D(w) — D(0) = 7, we obtain 



f (0) = j(g-2(0) ~ 9-2M) 



(52) 



If we set 7=1, the above equation reduces to equation 
(46) in LS. Equation (|52|) indicates a relation between 
the mass accretion rate and the second-order coefficient 
of density perturbation. The accretion rate is propor- 
tional to the difference of the second-order density per- 
turbation in the upstream and downstream flows and 
the proportionality constant depends on the velocity of 
the object as well as on the value of 7. Although equa- 
tion ([52j) gives us a physical insight of the accretion in 
terms of the coefficient <?_2, 7 and j3, we can not de- 
termine the accretion rate only based on this equation. 
Because we don't have a physical condition on the co- 
efficient of g_2 at = 0. Instead, we prescribe the 
accretion rate (i.e., /o(0)) from an accretion model. 

One should solve the two main equations (|48)) and 
(|49p numerically subject to the mentioned boundary 
conditions. We note that for the isothermal flow (i.e., 
7 = 1), our second-order equations (|48|) and (|49)) reduce 
to the equations (35) and (36) of LS. In order to solve 
the equations as a two point boundary value problem, 
the equations are integrated numerically from 9 = it 
to 6 = 0. Also, we need three initial values to start 
numerical integration from 9 = 0, two of which were 
introduced earlier, i.e. /o(tt) = fb(n) = 0. We have 
g_2 (tt) as the third boundary condition, but there is no 
constraint on its value. In order to obtain the physi- 
cal value of c/_2(7r), we follow the approach of LS, in 
which the mass accretion rate is prescribed via a phys- 
ical model to give us /o(0) = M. Then, we try to find 
g_2 (t) so that after solving the equations numerically, 
the prescribed value for /o at 9 = is obtained. The 
standard Bondi-Holye accretion model was used by LS 
in their model for determining the accretion rate. But 
we can not use the Bo ndi-Hoyle accretion rate, because 
our flow is isentropic. iFoglizzo fc Ruffertl (|1997|) did an 
extensive analysis to calculate appropriate relations for 
the accretion rate when the gas is isentropic. Their 
non-dimensional mass accretion rate is 



M{(3) = 2(i)^(^-i— )&(! + ( 



7-1 



)P 2 r- 




Fig. 5 The streamlines correspond to 7 = 5/3 (solid) and 
7 = 1 (dashed) for the cases with /3 = 0.5 and (7-2(1") = 1.12. 
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Fig. 6 The pressure contours correspond to 7 = 5/3 (solid) 
and 7=1 (dashed) for the cases with /3 = 0.5 and 3-2(71") = 
1.12. 



(53) 



Using equations (|53|) and (fTT|) , we obtain value of /o (0) 



for given j3 and 7 and it helps us to the calculate the 
proper value of g r _2(7r). 

The behaviour of the second-order coefficients /o (0) 
and <?_2 (0) are shown in Figures [TJ [21 and G3 For 
/3 = 0.5, we adopted three fiducial boundary values of 
9-2(71") (i.e., —0.88,0.12, 1.12) which allows us to make 
an easier comparison of our results for an iscntropic 
flow to the solutions of LS in the isothermal case. Each 
curve is labeled by the adiabatic index 7. Profiles of 
the second-order coefficients arc modified significantly 
when 7 varies from 1 to the larger values. Both fo(6) 
and 17-2(0) decrease for all values of g_ 2 (7r). The den- 
sity and the pressure contours and the streamlines are 
also shown in Figures 21 [5] and |B] corresponding to the 
two cases where 7 = 1 and 7 = 5/3. For a more physi- 
cal case, we used equation (|53p to obtain the accretion 
rate corresponding to the input parameters 7 = 1.4 and 
(3 = 0.5. Then, we have /o(0) = M = 1.25 and numer- 
ically it is found g_2(7r) = —0.15. Figure [7] shows the 
streamlines and the density contours for this case. 

Now, we are in a situation where we can derive the 
dynamical friction force for the isentropic problem. In 
order to do so, we calculate the rate of momentum in 
the z— direction that is transferred from the the object 
to the ambient medium. Considering a sphere of radius 
R centred on the object, the amount of the total mo- 
mentum transferred through the surface of this sphere 
is independent of its size as long as it lies outside the 
wake (LS). Thus, the total dynamical friction force con- 
sists of the kinetic and the static parts as follows: 



Fkin = — 27r / puRU z R sin I 
Jo 

and, 



F stat i C = -2?r / PR 1 cos sin 6d6. 
Jo 



(54) 



(55) 



Then, the non-dimensional dynamical friction force F 
becomes 



F 



gu T u z r sin9d9 — / p cos sin 



(56) 



Obviously, we have /i_ 2 (0) = 3-2(0) for an isothermal 
flow and so, the above relation for the friction force 
reduces to equation (57) of LS, i.e., 



- / ((l-/3 2 )sin6»cos 
/o 



0g_ 2 + 0(1 + cos 2 0)/J) del 



(58) 



LS showed analytically that equation (??) is intcgrablc. 
We also follow the same approach for evaluating the 
integral relation (??) for our isentropic flow. Like LS, 
the flow is irrotational and so the </>— component of the 
vorticity vanishes which leads to 



(3fo-(3 cot 9f' -P 2 sin 9 cos 9g'_ 2 +2(3 2 sin 2 6g- 2 = 



1-VD \ 
D 1 
(59) 



If we add the above equation to the equation (|48j) and 
multiplying by —(1/2) sin cos 0, we obtain 

sin 9 cos 0(/i-2-/? 2 .9-2)+/3 cos 2 9fo = - sin d cos 6H(D), 

(60) 

where H(D) is a function of D. Now, we can split 
equation (55) as 



f=~ r p/ d9 

Jo 

-A 

Jo 



sin cos 6{h- 2 - /3 2 .9- 2 ) + P cos 2 Bf ] d9. (61) 



Having equation (J5HJ), the second term of the above 
equation vanishes and we obtain 



F = pf (0) = PM, 



(62) 



where the unit of force is equal to 2ttpoc 2 R 2 . 

Now, we can substitute the corresponding expan- 
sions into equation (j5"6")l . We find that all terms with 
positive powers of r vanish upon integration, but the 
terms which arc independent of r survive. Therefore, 

F = - f [sin cos 9 (h _ 2 - f3 2 g _ 2 ) +/3(1 + cos 2 0) f a ]d6 
Jo 

(57) 



where M is the mass accretion rate in an iscntropic 
flow (i.e., equation (|53| ). Exactly the same relation is 
valid for the case of isothermal flows as demonstrated 
by LS. We can now study the effect of the adiabatic 
index 7 on the dynamical friction force in Figure [SI 
This figure displays the dynamical friction force as a 
function of the Mach number for different values of the 
adiabatic index. The dashed line shows the solution 
derived by LS for isothermal flows. This figure shows 
that the dynamical friction force decreases when 7 in- 
creases. Dependence of the dynamical friction force on 
the 7 becomes more significant when Mach number has 
a larger value. However, one should note that in our 
model our the accretion rate monotonically increases 

F'n the subsonic regime. But in the formula used by LS 
vhich is actually based on the Bondi-Holye model, the 
accretion rate monotonically decreases. This fact ac- 
counts for the divergence of all our curves from their 
result at intermediate Mach number. 
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Fig. 7 The streamlines and the density contours for a re- 
alistic case based on equation (|53[) for the accretion rate 
corresponding to 7 — 1.4 and j3 = 0.5. Here, we have 
/o(0) =M= 1.25 and g-2(n) = -0.15. 



In the present investigation, we considered a point-mass 
bject moving through an isentropic gas and the struc- 
jrc of the flow far from the object and the resulting 
ynamical friction force calculated scmi-analytically. In 
ict, our study is a direct generalization of the LS ap- 
roach for the isothermal gas to the isentropic case. 
Ising perturbation expansions for the density, pressure 
nd the stream function, we obtained the first-order 
nd the second-order differential equations describing 
le excited wakes when the object's velocity is sub- 
mic. We restricted our analysis to the steady state 
ase. In the subsonic regime, we proved analytically 
lat the dynamical friction force is MV, where M is 
lass accretion rate of the isentropic flow. Our results 
iow that a larger value of the adiabatic index corre- 
Donds to a reduction of the friction force for a given 
clocity. Our analytical solutions for the isentropic case 
snow that the wakes are less dense in comparison to the 
isothermal ones because the isentropic flow is less com- 
pressible than the isothermal case. This fact has also 
been pointed out by LS, but our study is the first to 
address this issue semi-analytically. 
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Fig. 8 The dimensionless friction force Fasa function of 
Mach number. Each curves is labeled by its value of 7. The 
dashed curve corresponds to the LS solution. 
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